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Abstract. - We study the properties of least time trajectories for particles moving on a two 
dimensional surface which consists of piecewise homogeneous regions. The particles are assumed 
to move with different constant speeds on different regions and on the boundary between regions. 
The speed of the particle is assumed to be highest when it moves along the edges formed by the 
boundary of two regions. We get an analogous behavior to Snell's Law of light refraction, but in 
a more generalized form. The model could be used for studying properties of animal and insect 
trails which tend to form predominantly along edges. The model predicts three types of behavior 
for the trajectories near a corner forming edge: fully edge following, partial edge following and 
complete avoidance of the edge, which are indeed observed in natural ant trails. 



Introduction. — The mathematical problem of find- 
ing the shortest or least time path has a long history. 
Heron of Alexandria (1 st century AD) was the first to real- 
ize that the law of reflection of light gives the shortest (or 
least time) path from a point to another on the same side 
of a dividing line which touches that line exactly once. In 
the seventeenth century Pierre de Fermat used the prin- 
ciple of least time to explain Snell's Law of refraction of 
light at the interface of two media in which the light trav- 
els with different speeds. The path that light ray takes 
in a medium with varying refractive index can be quite 
complex and is what is responsible for such phenomena as 
mirage and looming. 

Consider the refraction of light at the boundary of two 
media having different refractive indices. When one solves 
for the path by looking for the least time trajectory, the 
assumption is that the light either travels in one medium 
or the other but not in the interface of the two. But, it 
is possible to think of situations where the objects whose 
trajectories one is interested in might move along the in- 
terface as well, with a speed different from its value in the 
two media. The question to ask then is: How do the tra- 
jectories get modified with this additional freedom that 
the particle has? If we allow for particle motion along 
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the boundaries, the least time trajectories one gets are 
markedly different. The study of animal and insect trail 
patterns is one possible area where some of these trajec- 
tories might be relevant. 

Trail formation is a strategy that many organisms em- 
ploy for efficient execution of a collective task to be done. 
Trail formation in ants and termites are the most ubiqui- 
tous examples . Ants use well organized trails for trans- 
portation of food to nest as well as for shifting of colonies 
[l]. Another example of path optimization can be seen in 
the amoeba tube formation to the food source 2 . Trail 
formation is a fairly evolved collective endeavor that helps 
organisms to make the right choice (the most rewarding 
food source, for example), in the most economic manner 
(say, using the least 'costly' route from food source to 
nest) most of the times. It is thus interesting to model the 
shapes these trails tend to take. 

Ant trails typically tend to form along the shortest route 
from the nest to the destination, where for example there 
could be a food source 3 . But a notable exception to this 
behavior happens in the presence of structured guidelines 
or edges [2]. The ants then seem to go out of the way 
to accommodate appreciable amount of edges in its trail 
even when this would mean a longer route as compared 
to the straight line path from nest to food. The common 
example of this behavior being the trails formed inside 
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and outside the buildings along the wall edges, crevices, 
tile boundaries, table edges etc. It has been shown that 
in some species of ants, edge is used as a visual guide 
[5], whereas in other species they use it thigmotactically 
(by touch) [6]. One can model the behavior of ant trails 
by assigning different costs for travel on the surfaces and 
along the edges and trying to find the least cost trajectory. 
This is equivalent to the problem of finding the least time 
trajectories for particles moving with different velocities 
on different surfaces and along the edges, where speed of 
the particle is proportional to the inverse of the cost. 
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The model. — We shall study the particle motion on 
two dimensional surface which is piecewise homogeneous. 
That is, the surface is partitioned into separate regions in 
each one of which the particles move with a different speed. 
Additionally, we assume that the particle moves with a 
different speed if the motion is along the edges formed at 
the boundaries of two regions. The aim is to find the least 
time trajectories connecting two given points on such a 
surface. The problem for a most general distribution of 
regions is well defined but only numerically tractable. We 
shall look at couple of special and interesting cases where 
the trajectories are analytically solvable. 

Trajectories in the presence of a single straight 
boundary. — If there is only a single region present then 
the least time trajectory will be a straight line joining the 
starting to the destination point. Let us consider a sit- 
uation where two semi infinite regions (Si and S2) are 
separated by a straight boundary or edge (see fig. [I]). The 
starting point (N in fig. [T] ) is in Si and the destination 
point is in S2 (F in fig. 111. Let and Lp represent the 
perpendicular distance from the edge to N and F respec- 
tively. The horizontal distance between N and F is repre- 
sented by Lh- The particle velocities are v p when it moves 
along the boundary, and vi and V2 when it moves on the 
regions Si and S2 respectively. The quantities C{ = 1/vi 
can be thought of as the cost of traversing unit distance on 
the two regions (i = 1, 2) or the edge (i = p) and the least 
time problem can be equivalently posed now in terms of 
a least cost path. We shall take C p < C{ (i = 1, 2) as we 
assume that the particle has highest velocity when moving 
along edges. 

Let xi and X2 be the distances at which the particle en- 
ters and leave the edge respectively (see fig. [I]), measured 
from the point O in the figure. Let Oi be the angle the 
trajectory makes with the perpendicular at xi and O2 the 
angle the trajectory makes with the perpendicular at #2. 
The total cost associated with the trajectory, Cr, is 
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Here we are assuming that X2 is larger than x\. Minimiz- 
ing the cost Ct with respect to the parameter 61 and 62 



Fig. 1: Particle trajectory when the starting and the destina- 
tion points are separated by a straight boundary. Distances x\ 
and X2 are measured from the point O. 



we get, 
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This suggests that for fixed values of cost functions, entry 
angle in to the edge (0\) and the angle at which the par- 
ticle leaves the edge on to the second region (0 2 ) will be 
constants. 

The above analysis would break down if one ends up 
with a solution where X2 < #i, since in that case the cost 
would be negative for the portion along the edge (see eq. 
0). This would happen if 
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which can be obtained by noting that x\ = Livtan(#i) 
and X2 = Lh — LFtan(# 2 ) and using the equations ([2| 
and (|3|. In such a case there exists no solution in which 
the trajectory has the edge as its part. The correct cost 
function to be minimized in such a case would be, 
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which is obtained by putting X2 = X\ in eq. @. The 
solution that minimizes this cost function is given by 
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This gives the Snell's law for light refraction at the bound- 
ary of two surfaces. 

We see that the paths we have obtained is very similar 
to that of light ray refraction as the ray passes from one 
medium to another. The crucial difference in the present 
case being that the particle paths, unlike the light ray, can 
move along the interface of the two media (edge). The an- 
gle 0[ (6 2) never exceeds 61 (#2), which is like the upper 
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Fig. 2: Trajectories in the presence of corner forming bound- 
aries. The region on the convex side (region 2) has high cost 
of travel. (C2 > C\). The possible candidates for least energy 
trajectories are characterized by the two variables x v and xh- 



critical angle in the context of light refraction. But unlike 
in the case of light refraction where upper critical angle 
is present only for the ray that goes from a higher refrac- 
tive index medium to a lower refractive index one, in the 
present case upper critical angle is always present, as the 
speed travel along the edge is assumed larger than that 
for travel along the surfaces. 

Trajectories in the presence of a corner forming 
boundary. — We now study the path that the particle 
takes when edges formed between the two regions is not 
linear. We shall consider a right angled boundary formed 
when the surface is divided into two parts with one quad- 
rant being formed by type 1 region and the rest by type 2 
region as shown in fig. [2] We find the least cost trajectory 
when the source and the destination points lie on the two 
perpendicular boundary edges. We shall again take C p to 
be less than C2 and C\ . This problem has relevance to na- 
ture of ant trails formed on surfaces with corner forming 
guiding edges which we shall discuss in the next section. 

Taking the corner point (O in the fig. [2| as origin, Lh 
and L v are the distances to the source and destination 
measured from the origin. Since C p < C2 the minimum 
cost trajectory can have no part passing through region 2. 
The various possible trajectories can thus be parametrized 
using two points, Xh and x v . These are the point of depar- 
ture from edge on which source point lies on to region 1 
and point of entry into the perpendicular edge from region 
1. The minimum cost trajectory should have < Xh < Lh 
and < x v < L v . If C p and C\ represent the costs per 
unit length for traveling along the edges and along region 
1 respectively, the total cost incurred for a given path is 
given by, 

C b = C p (L h - x h ) + C P (L V - x v ) + \Jx 2 h + x 2 v Ci (7) 

where, < Xh < Lh and < x v < L v . 

Minimizing the cost C b with respect to the parameter 
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This tells us that if X = S 1 = y/2. there is an extremum 
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for all Xh = x v . The eigenvalues of the Hessian matrix 
evaluated at the points x h — x u are found to be and 
2* . The eigenvector corresponding to the eigenvalue 
is in the direction given by x v = Xh- This means that 
for X = y/2, there are a collection of paths all having the 
same costs and with x v = Xh- For other values of the 
ratio ^ no solution exists in the interior region defined 
by < Xh < Lh and < x v < L v . This implies that 
the minimum energy path in the Xh — x v plane should lie 
at the boundaries given by Xh = Lh , < x v < L v or 
x v = L v , < Xh < Lh or Xh = x v = 0. Note that if 
either x v or Xh is zero, the trajectory corresponds to the 
one fully along the edge. 

In order to find the minima when X =fi y/2 we separately 
minimize the cost on the boundary regions given above. 
The cost incurred when x v = L v is given by, 



C' b = C p (L h - x h ) + + LI & . (9) 

Minimizing the above cost function with respect to Xh we 
get, 
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where C f min is the corresponding minimized energy. Sim- 
ilarly, substituting Xh = Lh in Eq. Q and minimizing 
the cost function with respect to x v one gets the possible 
paths which minimizes the cost function given by, 
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where C^ in is the corresponding minimized energy. The 
cost incurred when the particle moves completely on the 
edge (x h = x v = 0) is C e min = C P (L V + L h ). 

The trajectory will deviate from the edge if one of the 
minimum costs computed above, C^ in or C^ in , is less 
than the cost incurred while traversing the trail fully along 
the edge, C^ in . The condition for the trail to deviate 
from the edge is thus given by, ^ < y/2. Further, if 
we assume that Lh > L v (the other case can be treated 
similarly), the solution for which Xh = Lh is not valid 

because x Vl 



Lh/ y Cf/Cp — 1 becomes larger than L v . 
This definitely cannot be a minima. This spurious minima 
turns up because the cost function that we are using is 
not appropriate when x v > L v or Xh > Lh- Thus for 
Lh > L v , the solution will have x v = L v and < Xh < Lh- 
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partly on the surface 



Purely edge following 
behavior 




Fig. 3: The phase diagram showing various kinds of trajectories 
formed in the presence of corner forming edges assuming Lh > 
L v 



For values of X less than X = y/l + L^/L^ the value for 
Xh that minimizes the energy will become larger than Lh 
and thus unphysical. This implies that for those values 
of Ci/Cp, the minimum energy path would be given by 
x h L h and x v L v . 

The results we have derived above is shown in the phase 
diagram (see fig. |3j). The phase diagram shown is for 
the case when Lh > L v . For values of X above y/2, the 
minimum energy path is the one where the trail is along 
the edge completely. For values of X lying between y/2 and 
y/l + L^/L^, the minimum energy path is partly along 
the horizontal edge and the rest on region 1. For values 
of X less than y/l + L%/L^, the most economic path is 
a straight line connecting source to destination. In the 
limit Lh ^> L v , the range of X values over which one can 
observe a path that is partially on the edge and partially 
on the surface is maximum. In the limit of Lh = L v , the 
path is either completely along the edge or completely on 
the surface depending on whether the value of X is larger 
that y/2 or not. 

It is possible to generalize the above analysis to the case 
where the starting and the destination points do not lie on 
the edges themselves. In this case the trajectories which 
has edges as its parts will then have to be characterized 
with more number of parameters and one has to compare 
the cost of such paths with the cost of the straight line path 
between the two points. We shall discuss below certain 
observations made on natural ant trails where the model 
and the particular case we have looked at can qualitatively 
explain the features seen. 

Qualitative comparison with naturally observed 
trails. — We have made some preliminary observa- 
tions on the edge following feature of ant trails (see figs. 
|4(a)|4(b)|4(c) and |4(d)| ) and they are qualitatively found 
to exhibit the types of behavior that has been seen in the 
analysis above. The ant trails have been photographed 
and then traced. Even though the locations of the ants 
come with a small error (zblcm) due to visual copying of 
the ant locations from the photographs on to the figures, 
the overall feature of the trails are accurate. Some of the 
trails given in the figures were formed naturally, but in 



couple of the cases (fig. 4(a) and fig. |4(d)[ ), we intro- 
duced the food sources to let the trail form in the regions 
of interest. 

The most common situation that one encounters is when 
there is a food source on the floor, close to a straight edge 
at the intersection of a wall and the floor. The ant trails 
will typically take off at an angle from the edge and ap- 
proach the food in a straight line as seen in fig. |4(a)l The 
situation here is much like the case that we have initially 
analyzed (see fig. [T]) above with the starting point lying 
on the boundary. The most economic path that the ants 
choose is not the shortest route but a route that partially 
uses the edge (boundary) for navigation. Let us now look 
at a few trail patters in the presence of corner forming non 
straight edges. In fig. 4(b) we have shown an ant trail that 
was seen to avoid the edges completely and form a straight 
line path from the nest to destination directly. The des- 
tination in the present case was the corner of the ledge, 
beyond which the trail was along the ledge in a direction 
perpendicular (into the paper in the figure) to the surface 
shown. Fig. 4(c) | shows the case where the trail stuck to 
the edge completely while navigating the corner. This is 
one of the phases (when X > y/2) that is expected in the 
analysis. The trail shown in fig. |4(d)| is typical of a situ- 
ation where a long wall meets a short one leading to two 
perpendicular edges meeting with Lh ^> L v (in the present 
case L v = 9cm and Lh ~ 50cm). The ant trail is seen 
to have formed partially on the edge andjpartially on the 
surface, conforming to the phase in fig. |3]when X < y/2 
and X > 1. Occasionally, one also comes across trails 
where x v is not equal to L v or (with Lh > L v ), which 
contradicts our analysis. This could be due to the fact 
that the local rules that ants employ to find the minimum 
cost trail does not always lead to the path with least cost. 

Ideally, an independent calculation of X should be done 
to make a comparison between the model and the experi- 
mental data. That would require experiments to be done 
with complete control over the type of cues that ants could 
use in the trail construction and maintenance. There are 
indications of such results in experiments that were done 
in the context of nocturnal orientation of black carpenter 
ants by Klotz et al. [7]. In the arena test conducted in 
this work, it was seen that shielding visual cues progres- 
sively lead to more and more edge following thigmotactic 
behavior. Experimental verification of the phase diagram 
we are proposing can be done using similar experimental 
settings. Another instance where a Snell's law like behav- 
ior has been seen is in the context of light avoiding amoeba 
path formation towards food source [2]. 

Summary. — We have started with a very simple 
model, reminiscent of Hamilton's principle for light rays, 
to study the nature of particle trajectories in a surface 
formed by different regions on which the particles move 
with different speeds. The particles are assumed to move 
along the edges formed by the intersecting surfaces with a 
different speed. This leads to a generalized form of Snell's 
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Fig. 4: Different types of ant trails observed near wall edges and corners. The shaded regions in the figures correspond to 
surfaces coming out of the plane of the paper, (a) The trail formed along an edge cuts at an angle (« 24°) to the edge to access 
the food lying on the floor, (b) The ant trail completely avoids the edges and forms directly from the nest to the destination 
(Lh = 13cm, L v = 9cm). (c) The case where the trail was seen to follow the edge completely, (d) Trail is formed partially along 
the edge and partially on the surface (L v = 9cm, Lh = 50cm). 



Law in the context of these particle paths. In the case 
where the boundaries are not straight, the analysis shows 
that non trivial trajectories can be formed: from com- 
pletely edge following to partial edge following and com- 
plete avoidance of the edges. Our observations of real ant 
trails confirms the existence of these various phases. The 
comparisons we have done are qualitative but quantita- 
tive study of the phase diagram we find is possible. In 
the context of ant trails, we have not discussed about the 
methods by which ants could be zeroing in on the optimal 
trail. This has to be a set of local rules that ants employ 
to arrive at better routes 8,9. How these local rules are 
affected in the presence of linear or piece wise linear guide- 
lines will be an interesting question to address. The ideas 
presented here can be of use in autonomous vehicle design 
which could use edges for navigation and for optimizing 
the functions of chemical trail laying ant robots 10 . That 



one could be seeing glimpses of something as fundamental 
as Snell's law, in fact a more generalized version of it, in 
something as ubiquitous as ant trails is both surprising 
and interesting. 
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